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ABSTRACT 

We use recent HST colour-magnitude diagrams of the resolved stellar populations 
of a sample of local dSph galaxies (Carina, Leol, LeoII, Ursa Minor and Draco) to 
infer the star formation histories of these systems, SFR{t). Applying a new variational 
calculus maximum likelihood method which includes a full Bayesian analysis and al- 
lows a non-parametric estimate of the function one is solving for, we infer the star 
formation histories of the systems studied. This method has the advantage of yielding 
an objective answer, as one need not assume a priori the form of the function one is 
trying to recover. The results are checked independently using Saha's W statistic. The 
total luminosities of the systems are used to normalize the results into physical units 
and derive SN type II rates. We derive the luminosity weighted mean star formation 
history of this sample of galaxies. 

Key words: methods: statistical - stars: formation - galaxies: evolution - Local 
Group 



1 INTRODUCTION 

■ The local dwarf spheroidal galaxies form a sample of small 
'■ galaxies which, due to their relatively nearby locations and 
. close association with the Milky Way could in principle fur- 
' nish crucial observational and theoretical information on 
a range of astrophysical phenomena. As these systems are 
eventually disrupted and incorporated into the Milky Way 
they illustrate locally one of the mechanisms thought to be 
responsible for the build up of large galaxies. Thus com- 
paring their stars and those now found in our galaxy we 
can obtain a first estimate of the relevance of late mergers 
(Unavane et al. 1996). Dynamical studies of their stars have 
yielded valuable constraints on the nature and structure of 
dark matter halos at the smallest scales (Lin & Faber 1983, 
Gerhard & Spergel 1992). Their orbits probe the galactic 
halo in a range of distances not sampled by any other ob- 
jects and can thus be used to study the outer galactic halo. 
Additionally, their small sizes make them in principle the 
simplest galactic systems, where key processes such as star 
formation and gas flows can be studied under relatively well 
defined conditions. 

However, this situation is complicated by the fact that 
we not only lack a theoretical understanding of these sys- 
tems, but also an observational record of their evolution; 



only a present day snapshot of their physical parameters is 
available, as is the case with most galactic systems. Whilst 
high redshift observations have recently opened up new areas 
of research as they begin to yield an statistical description 
of the evolution of bright galaxies, such an approach is likely 
to remain out of reach for these small systems for some time. 
Fortunately, their neighboring locations allow the study of 
their individual stars, which ofi^ers the possibility of directly 
probing their evolutionary histories by inferring star forma- 
tion rates as a function of time, SFR{tYs. 

The recent availability of detailed colour-magnitude di- 
agrams for several nearby systems has prompted the de- 
velopment and allowed the application of careful statistical 
methods aimed at reconstructing the star formation histo- 
ries of these objects (e.g. Chiosi et al. (1989), Aparicio et 
al. (1990) and Mould et al. (1997) using Magellanic and lo- 
cal clusters, and Mighell & Butcher (1992), Smecker-Hane et 
al. (1994), Tolstoy & Saha (1996), Aparicio & Gallart (1995) 
and Mighell (1997) using local dSph's). Although much has 
been learnt of the complex SFR{tYs of these systems, ex- 
isting studies have lacked two major ingredients: a homo- 
geneous set of observations including several of the dSph 
galaxies does not exist, and different data sets are gener- 
ally analyzed using different techniques. These two points 
make comparisons between the derived SFR{tYs at best 
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risky. A further difficulty lies in the fact that the available 
rigorous statistical studies approach the problem paramet- 
rically, which is something one should try to avoid when 
the actual structure of the function one is trying to recover 
can be crucial, as is the case when the underlying physics 
is unknown. An example of this last point is the case of 
the Carina dwarf. Hurley-Keller et al. (1998) solve for the 
best fitting three discrete bursts solution to the SFR{t) and 
conclude that star formation has proceeded spasmodically, 
whilst MighcU (1997) uses a non parametric star count ap- 
proach, albeit not a fully consistent statistical method, to 
obtain a more gradual solution for Carina's SFR{t). 

In this paper wc have attempted to improve on the de- 
termination of the star formation histories of local dSph sys- 
tems by addressing the two points mentioned above. Wo 
use recent HST observations of the resolved populations of 
a sample of dSph galaxies (Carina, Leol, LeoII, Ursa Mi- 
nor and Draco) uniformly taken and reduced, to recover the 
SFR{t) of each, applying a new non-parametric maximum 
likelihood method. This allows meaningful comparisons to 
be made, as any systematics, at any level, will affect all our 
galaxies equally. 

The outline of our paper is as follows: in section 2 wc 
discuss the observations, in section 3 we include a brief out- 
line of our method, which was introduced in our paper I 
(Hernandez et al. 1999). The results axe presented in sec- 
tion 4, and in section 5 we summarize our results. 



2 THE OBSERVATIONS 

The main requirements of our observations were that they 
should comprise a homogeneous sample of local dSph galax- 
ies, mostly in terms of the data reduction. Only such an 
internally consistent data set allows robust comparisons 
between different galaxies, once uniform data reduction 
and analysis methods are adopted. We extracted available 
archive HST data for the Carina, Leol, LeoII, Ursa Minor 
and Draco galaxies, and used standard data reduction meth- 
ods and standard HST calibration numbers throughout the 
sample (e.g. Elson et al. 1996). 

The currently available data cover only small (and vari- 
able) sections of the total extent of these systems. This fact 
clearly limits the inferences which can be drawn to the small 
observed fractions, the star formation histories of these re- 
gions might not by representative of the average for a whole 
galaxy. While this limitation introduces an extra uncertainty 
to our results, it highlights the interesting possibility of 
studying spatial variations in the evolutionary histories of 
dSph galaxies, if comprehensive HST studies where under- 
taken. In the above sense, our results for the different galax- 
ies refer in the strictest sense only to the fractions covered 
by the observed fields. 

As we did not require the HR diagrams to extend much 
fainter than the oldest turnoff points (My « 24 — 25), or 
to be complete into the faintest limits (the faintest stars 
were in fact excluded from the analysis) the data reduction 
was straightforward. Our resulting CMDs do not show any 
systematic difference from comparable published ones for 
the galaxies we study. The technical details of the images 
used appear in the appendix. 



3 THE METHOD 

In this section we give a summary description of our HR di- 
agram inversion method, which was described extensively in 
our paper I. In contrast with other statistical methods, we do 
not need to construct synthetic colour magnitude diagrams 
(CMD) for each of the possible star formation histories be- 
ing considered. Rather we use a direct approach which solves 
for the best SFR{t) compatible with the stellar evolutionary 
models assumed and the observations used. 

The evolutionary model consists of an isochrone library, 
and an IMF. Our results are largely insensitive to the details 
of the latter, for which we use: 
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The above fit was derived by Kroupa et al. (1993) for a 
large sample towards both galactic poles and all the solar 
neighborhood. 

As the weak metallicity dispersions measured in the 
galaxies we are studying (around 0.3 dex) are comparable to 
the errors in the metallicity determinations themselves, we 
have not attempted to introduce any enrichment histories 
for any of our galaxies. In fact, as small internal metallic- 
ity spreads present in these galaxies would introduce only 
small differential age offsets in our inferred SFR{t) (see Ta- 
ble I), we shall in all cases use single metallicity isochrone 
sets. Once a metallicity has been selected, wc use the latest 
Padova isoclirones (Fagotto et al. 1994, Girardi et al. 1996) 
together with a detailed constant phase interpolation scheme 
using only stars at constant evolutionary phase, to construct 
an isochrone library having a chosen temporal resolution. 

In this case we implement the method with a resolution 
of 0.15 Gyr, sufficient for our present problem. It is one 
of the advantages of the method that this resolution can 
be increased arbitrarily (up to the stellar model resolution) 
with computation times scaling only linearly with it. 

Our only other inputs are the positions of, say n ob- 
served stars in the HR diagram, each having a colour and 
luminosity, c; and k. Starting from a full likelihood model, 
we first construct the probability that the n observed stars 
resulted from a certain SFR{t). This will be given by: 



=n(/ 

i=l ^"''o 



SFR(t)Gi{t)dt 



(2) 



where 



Gi{t) = 



p(U:t) 



2-Ka(li) 



exp 



In the above expression p{li; t) is the density of points 
along the isochrone of age t, around the luminosity of star 
i, and is determined by the assumed IMF together with the 
duration of the differential phase around the luminosity of 
star i. to and ti are a maximum and a minimum time needed 
to be considered, for example and 15 Gyr. a{li) is the am- 
plitude of the observational errors in the colour of the stars, 
which are a function of the luminosity of the stars. This 
function is supplied by the particular observational sample 
one is analyzing. Finally, C{li\t) is the colour the observed 
star would actually have if it had formed at time t. We shall 
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refer to Gi{t) as the likelihood matrix, since each clement 
represents the probability that a given star, i, was actually 
formed at time t. Since the colour of a star having a given 
luminosity and age can sometimes be multi-valued function, 
in practice wc check along a given isochrone, to find all possi- 
ble masses a given observed star might have as a function of 
time, and add all contributions (mostly 1, sometimes 2 and 
occasionally 3) in the same Gi{t). In this construction we 
are only considering observational errors in the colour, and 
not in the luminosity of the stars. The generalization to a 
two dimensional error ellipsoid is trivial, however the obser- 
vational errors in colour dominate the problem to the extent 
of making this refinement unnecessary. Although the ampli- 
tude of luminosity errors is only a factor of ~ 2 smaller than 
colour errors, as can be inferred from the fact that CMD 
diagrams typically display a range of luminosities 5 times 
larger than in colour, in discriminating between isochrones, 
errors in colour are « 10 times as important as errors in 
luminosity. The absence of a colour dependence from p{li; i) 
is a direct consequence of having neglected errors in the lu- 
minosity of the staxs. A star of a given observed luminosity 
and assumed age will thus have a colour determined by the 
isochrones used. 

Equation (2) is essentially the extension from the case of 
a discrete SFR{t) used by Tolstoy & Saha (1996), to the case 
of a continuous function (continuous in time, but obviously 
discrete with respect to the stars) in the construction of the 
likelihood. The challenge now is to find the optimum SFR{t) 
without evaluating equation (2) i.e. without introducing a 
fixed set of test SFR{t) cases from which one is selected. 

The condition that jC{SFR) has an extremal can be 
written as 

5C{SFR) = 0, 

and a variational calculus treatment of the problem applied. 
Firstly, we develop the product over i using the chain rule 
for the vaxiational derivative, and divide the resulting sum 
by C to obtain: 



E 



= 



(3) 
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Introducing the new variable Y{t) defined as: 

Y{t) = j ^JSFR{t)dt =^ SFR{t) = (^^^^^ 

and introducing the above expression into equation (3) 
we can develop the Euler equation to yield. 
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where 



= / SFR{t)Gi{t)dt 

J to 

This in effect has transformed the problem from one of 
searching for a function which mgiximizes a product of in- 
tegrals (equation 2) to one of solving an integro-difFerential 
equation (equation 4). We solve this equation iteratively. 



Table 1. Age offset in Gyr between simulated and recovered pop- 
ulations as a function of input age (in Gyr) and metallicity (in 
solar units), for a ± 0.2 dex metallicity mismatch 



[Fe/H] 


-1.4 


-1.6 


-1.8 


-2.0 


-2.2 


2 Gyr 


-t-0.6 


+0.4 
-0.3 


+0.3 
-0.1 


+0.4 
-0.2 


-0.2 


4 Gyr 


-1-0.8 


+0.5 
-0.6 


+0.4 
-0.4 


+0.6 
-0.2 


-0.4 


6 Gyr 


+1.6 


+0.5 
-0.9 


+0.7 
-0.2 


+0.9 
-0.5 


-0.6 


8 Gyr 


+1.8 


+0.4 
-1.3 


+0.5 
-0.8 


+1.4 
-0.6 


-0.9 


10 Gyr 


+2.5 


+0.5 
-2.0 


+0.7 
-0.5 


+1.3 
-0.7 


-1.5 



with the boundary condition SFR(15)=0. Details of the nu- 
merical procedure required to ensure convergence to the 
meiximum likelihood SFR(t) can be found in our paper I, 
where the method is tested extensively using synthetic HR 
diagrams. The main advantages of our method over other 
maximum likelihood schemes are the totally non parametric 
approach the variational calculus treatment allows, and the 
efficient computational procedure, where no time consuming 
repeated comparisons between synthetic and observational 
CMD are necessary, as the optimal SFR{t) is solved for di- 
rectly. 

The lower main sequence region of the CMD diagram is 
totally degenerate with age, and contains the lower bright- 
ness stars, where the errors are larger. We have seen from 
using synthetic HR diagrams that excluding this region pro- 
duces a faster and more accurate convergence of the method, 
and have in analyzing real galaxies excluded stars of magni- 
tudes fainter than My ^ +5. This lEist cut together with the 
fact that our isochrones only extend out to the tip of the red 
giant branch (to go further would necessitate combining re- 
sults from different physical models, which we preferred not 
to do) leaves us with a mass range which actually varies as a 
function of time. To include also the fraction of the SFR{t) 
outside this region, we apply a minor correction factor to 
the result of equation (4), which accounts for the fraction 
of mass outside the sampled range, as a function of time, as 
given by the IMF used. 

Before presenting the star formation histories which re- 
sult from applying our method to the colour-magnitude di- 
agrams of the galaxies sampled, we include a summary of 
the systematic errors associated with the theoretical inputs 
of the method, which are more extensively discussed in our 
paper I. 

The IMF convolved with the duration of the differential 
evolutionary phase enters the calculation of the likelihood 
matrix in determining the density of points around the lumi- 
nosity of each observed star, for each of the isochrones con- 
sidered. As the main sequence and the giant branch regions 
of the CMD are degenerate with age (for a single metallicity 
population), it is the region containing the turn-off points 
of the sampled population that drives the solution of the 
problem. As a consequence of the above, the details of the 
IMF used are largely unimportant, it is basically the main 
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sequence lifetime of a star that the solution is sensitive to. 
This was shown explicitly in our paper I through the use 
of synthetic HR diagrams, where the IMF affected only the 
amplitude of the recovered SFR{tYs, which were normal- 
ized through the total number of stars in the HR diagrams. 
As in this case we are normalizing the inferred SFR{tys 
through the total luminosities of the galaxies being studied, 
changing the IMF within any reasonable limits leaves the re- 
sults unaffected. The effect of any blended binaries is equally 
unimportant, as the broadening of the main sequence occurs 
in the degenerate region and is in any case much smaller 
than the broadening produced by the observational colour 
spread. 

The resolution of the method varies as a function of 
time and of the observational errors present in the CMD be- 
ing analyzed, in ways that were studied in our paper I. The 
observational errors tend to smear the time structure in the 
SFR{t) always towards older ages, i.e. a burst of ago t will 
be recovered as an episode of duration At ending at time 
t. This At varies with age, becoming significant (> IGyr) 
only for populations older than around 10 Gyr, for the level 
of observational errors present in our HST colour-magnitude 
diagrams. Younger populations are less affected, and the for- 
mal resolution at which the method was implemented of 0.15 
Gyr is representative of our results for ages younger than 6 
Gyr. 

As explained in the previous section, our isochrones end 
at the tij) of the red giant branch, which means that more ad- 
vanced evolutionary phases can not be incorporated into the 
analysis. Fortunately, these later phases occupy regions of 
the CMD diagram distinct from those containing the phases 
we account for. Therefore, we can simply remove any red 
clump and horizontal branch stars from the analysis, leav- 
ing the structure of the studied regions unaifected. These 
later phases form only a minority component, containing 
little extra information, and do not affect our inferences. In 
the same way our results arc not affected by the presence of 
some contaminating field stars. Provided they do not fall on 
the region of the CMD containing the turn off points of the 
underlying stars, they are simply removed from the analysis. 

We are however sensitive to the assumed metallicity of 
the stellar populations being treated. In our paper I we pre- 
sented a few examples of how the method reacts when in- 
verting a synthetic HR diagram produced with isochrones 
different to the ones used in the inference procedure. If we 
construct the likelihood matrix using isochrones which arc 
very different (about 1 dex) from the ones used to produce 
the HR diagram, the iterative method used in solving equa- 
tion (4) becomes unstable and tends to divergent solutions. 
This property can be used to deduce large incompatibilities 
between the stars and the template isochrones against which 
they are being compared. Small metallicity offsets are much 
harder to detect, and produce distorted results. The degree 
of distortion varies with age, shape of the overall SFR{t) and 
the observational errors present in highly non-linear fashion. 

To give some indication of these distortions we present 
Table I. We produced a synthetic HR diagram from a sin- 
gle Gaussian burst input SFR{t) having a duration of 1 
Gyr, of a given metallicity, and applied the method using a 
metallicity 0.2 dex lower than the one used to generate the 
stars. This was repeated for a range of metallicities and ages, 
for both positive and negative metallicity mismatches. The 



top row of Table 1 shows the input metallicity (in dex) , the 
first column shows the input burst age, and the other entries 
show the age offset between the recovered and input SFR{t), 
all ages arc given in Gyr. A "-f" sign denotes the inferred 
population was older than the input one, where a lower 
metallicity was used in the inference procedure. Similarly, a 
"-" indicates that the inferred population was younger than 
the input one, where a higher metallicity was used to invert 
the simulated CMD. 

The metallicities shown cover the range present in the 
galaxies being studied. Populations much younger than 2 
Gyr are not present in our galaxies, and those older than 10 
Gyr arc distorted by the observational errors to the point of 
eliminating much of the time resolution of the inversion in 
this region. This table can be used to estimate the effects of 
changing the assumed metallicities, or of introducing tem- 
poral gradients, within the small observationally restricted 
range. The well known age-metallicity degeneracy is appar- 
ent. This might affect our results even if the mean metallic- 
ity is well known, as temporal variations in the metallicities 
(which must exist at some level) are not considered by the 
method. Using an independent test on our results, we find 
these effects to be minor in most of the cases we study, as ob- 
servational measurements of the metallicity of these systems 
suggest. 

3.1 Testing the results 

Once the IMF, metallicity and observational parameters are 
assumed for a given galaxy, the positions of the observed 
stars in the CMD are used to construct the likelihood ma- 
trix Gi{t), which is the only input given to the inversion 
method. In our Paper I we tested this method using syn- 
thetic CMDs produced from known SFR{t)' with which 
we could assess the accuracy of the result of the inversion 
procedure. In working with real data, we require the intro- 
duction of an independent method of comparing our final 
result to the starting CMD, in order to check that the an- 
swer our inversion procedure gives is a good answer. From 
our paper I we know that when the stars being used in the in- 
version procedure were indeed produced from the isochrones 
and metallicity used to construct the likelihood matrix, the 
inversion method gives accurate results. The introduction 
of an independent comparison between our answer and the 
data is hence a way of checking the accuracy of the input 
physics used in the inversion procedure, i.e. the IMF, metal- 
licity and observational parameters. 

The most common procedure of comparing a certain 
SFR{t) with an observed CMD is to use the SFR{t) to 
generate a synthetic CMD, and compare this to the ob- 
servations using a statistical test to determine the degree 
of similarity between the two. For example, Apaxicio et al. 
(1997) manage to recover simultaneously the distance, en- 
richment history and SFR{t) of the local dwarf LGS 3 by 
constructing synthetic CMD's from a SFR{t) taken para- 
metrically as a series of contiguous bursts, and finding the 
amplitudes of each burst that give the maximum statistical 
similarity with the data, in terms of a counts in cells maxi- 
mum likelihood. In that case, they solve for the amplitudes 
of the bursts that give a total synthetic CMD most closely 
resembling the observational data set. Their synthetic CMD 
is a linear sum of the partial CMD's produced from a single 
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Figure 1. 

Left: Observational HR diagram for Carina. Right: Inferred SFR{t) for the central values of the observational parameters, solid line. 
The dashed curves represent the error envelope as defined by the quoted uncertainties in the foreground extinction and distance 
modulus. 



realization for each burst. This has the advantage of allowing 
a large parameter space to be considered, as the synthetic 
CMD's are constructed trivially from a fixed single statis- 
tical realization of each burst. The disadvantage however is 
that one is not comparing the SFR{t) with the data, but 
rather a particular realization of the SFR{t) with the data. 
The distinction becomes arbitrary when large numbers of 
stars are found in all regions of the CMD, which is generally 
not the case. Following a Bayesian approach, we prefer to 
adopt the W statistic presented by Saha (1998), essentially 

-pr K + aQ! 

i = l 

where B is the number of cells into which the CMD is 
split, and rrii and Si are the numbers of points two distri- 
butions being compared have in each cell. This asks for the 
probability that two distinct data sets are random realiza- 
tions of the same underling distribution. In implementing 
this test we first produce a large number (~ 500) of random 
realizations of our best answer SFR{t), and compute the W 
statistic between pairs in this sample of CMD's. This gives 
a distribution which is used to determine a range of values 
of W which are expected to arise in random realizations of 
the SFR{t) being tested. Next the W statistic is computed 
between the observed data set, and a new large number of 
random realizations of SFR{t), this gives a new distribu- 
tion of W which can be objectively compared to the one 
arising from the model-model comparison to assess whether 
both data and modeled CMD's are compatible with a unique 
underling distribution. Both distributions of W were char- 
acterized in terms of a mean value and a la amplitude. This 



final check of our answer is in fact the slowest part of the 
procedure, but necessary to obtain an independent check on 
the answer of our inversion method. In other terms, we are 
checking that our best inferred maximum-likelihood solution 
is also a good fit. The value of B used was ~ 6400. 



4 THE GALAXIES 
4.1 Carina 

The first galaxy we study is the Carina dwarf, one of the 
first dSph's to be observed in terms of resolved stellar pop- 
ulation, and the one for which the most studies inferring 
the SFR{t) from the CMD have been published. This gives 
us the opportunity of comparing our results with previous 
studies. Our CMD diagram initially contained 2550 stars. 
After removing contamination, stars beyond the RGB and 
the lower degenerate region, we are left with 980 stars. The 
full observational CMD is shown in the left panel of Figure 
(1). 

Using the numbers published in the recent review by 
Mateo (1998) we took as the central values of our observa- 
tional parameters for Carina [Fe/H] — —2.0 ± 0.2, Mv = 
-9.3, E{B -V) = 0.04 ± 0.02 and (m - M)o = 20.03 ± 0.09 
(Smecker-Hane et al. 1994, Mighel 1997, Hurley-Keller et 
al. 1998 and Mateo et al. 1998). The metallicity fixes the 
isochrones and the integrated magnitude the normalization, 
we use the distance modulus and the reddening correction 
to fix the observations in the CMD. At this point we apply 
our method to invert the observational CMD and recover the 
underlying SFR{t), this is shown in the right panel of Figure 
(1) by the solid curve. The dotted curves represent the upper 
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Figure 2. 

Left: One synthetic HR diagram for Carina produced using the central inferred SFR{t) of Figure (1). Right: Distribution of values of 
W arising from 500 model-model comparisons, unshaded histogram. This establishes the range of values of W which result from random 
realizations of the inferred SFR{t). Also shown is the distribution of values of W arising from 250 data-model comparisons, showing 
the hypothesis that the data result from the same underlying function as the models to be acceptable at more than 1 sigma level. 



and lower envelopes to a series of alternative reconstructions 
of the SFR{t) produced by changing the assumed values of 
the distance modulus and the reddening corrections, within 
their respective error ranges. The internal metallicity spread 
present in this galaxy is very low, as can be seen from the 
narrow RGB, and is quoted as < 0.1 dex in the review by 
Mateo (1998). In this way, our error margins for this galaxy 
due to metallicity uncertainties are not much larger than 
what is shown by the dotted lines, see Table(l). 

In Figure (2) we illustrate the procedure of checking the 
inferred SFR{t), the left panel gives one random realization 
of the central inferred SFR(t), which is seen to resemble 
the data for Carina rather well. The right panel of Figure 
(2) shows the implementation of the W test. The solid his- 
togram gives the distribution of values of W which arise 
from 500 model-model comparisons, and gives the variability 
arising from the different random realizations of the central 
SFR{t), for the number of stars present in our observations. 
The dashed histogram gives the distribution of values of W 
which result from 250 data-model comparisons. Only < 32 
% of the random realizations of our central SFR{t) would 
give distributions of W (when compared against all other re- 
alizations) having a mean value further removed from that 
of the 500 model-model distribution than the data. In this 
sense, we can accept the hypothesis that both the data and 
our 500 random realizations of the central SFR{t) for Ca- 
rina come from the same underlying generating function at 
a 1 0" level. For the remaining galaxies we shall give only the 
results of the W test in terms of the mean and la amplitude 
of the model-model and the data-model distributions. 

Our result shows an interesting SFR{t) for this galaxy. 



very little star formation at early times, until around 10 
Gyr ago, when over a period of 3 Gyr an intermediate pop- 
ulation was formed. The SFR{t) then decreased markedly, 
before entering a more recent and extended period of star 
formation which ended 2 Gyr ago. The very low levels found 
throughout for the star formation rate of 50 — IOOMq /Myr 
are representative of what we find in all our galaxies, and 
should provide clues as to the physical processes driving the 
star formation activity in these systems. 

The existence of some RR Lyr stars in this galaxy (e.g. 
Saha et al. 1986, Mateo et al. 1995 or Kuhn et al. 1996) signal 
the presence of a very old population at some level, although 
recent estimates of the SFR{t) in Carina coincide in that the 
amplitude of this old component is minor (Hurley-Keller & 
Mateo 1998), and appears to have blended completely into 
the MS of our CMD, or to be found preferentially in a region 
of the galaxy not sampled by our observations. One of the 
most general claims about the SFR{t) of Carina has been 
the extreme "bursting" character of it (e.g. Smecker-Hane 
1994, Hurley-Keller & Mateo 1998) we note however, that 
all these studies have assumed a priori an extremely dis- 
crete form for the SFR{t) of this galaxy, and then solved 
for the best such function. In contrast, Mighell (1997) using 
a non-parametric approach, finds a much more continuous 
solution, basically consistent with what we obtain. The res- 
olution of our method in the region of ages 2 — 8Gyr is 
sufficient to exclude the possibility of any total cessation of 
the star forming activity in this region lasting more than 
0.5Gyr, as was shown in our Paper I, where synthetic HR 
diagrams produced from bursting SFR{tYs were correctly 
inverted. We conclude that although the star formation his- 
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Figure 3. 

Left: Observational HR diagram for Ursa Minor. Right: Inferred SFR(t) for the central values of the observational parameters, solid 
line. The dashed curves represent the error envelope as defined by the quoted uncertainties in the foreground extinction and distance 
modulus. 



tory of this galaxy is clearly bimodal, it is not a series of 
discrete bursts lasting < IGyr. However, we can not ex- 
clude the possibility that sampling a much larger fraction of 
this galaxy could yield somewhat different results. Analyz- 
ing local variations of the SFR{t) within these galaxies is 
an interesting project which will be treated in other papers, 
using data covering greater portions of the galaxies. As with 
all our galaxies, given the small number of stars available, 
the resolution of our solution is limited. Episodes of very 
short duration or low level, resulting in very few stars will 
be totally missed. We are recovering the SFR{t) responsible 
for producing the greater fraction of the observed stars. 



4.2 Ursa Minor 

The case of Ursa Minor seems to be the simplest of the ones 
we study, and is actually the only one of our galaxies which 
agrees with the once common expectation of dSph systems 
being simply old and metal poor. Our observational CMD 
is shown in the left panel of Figure (3), and is made up of 
1232 stars. After removing those stars incompatible with the 
phases included in our isochrones, together with the fraction 
beneath My = 6 we are left with 334 stars, which we used 
in the inversion procedure. 

Using the values given in Mateo's (1998) review, we take 
[Fe/H] = -2.2 ± 0.1, Mv = -8.9, E{B ~ V) = 0.03 ± 0.02 
and (m - M)o = 19.11 ± 0.1 (Nemec et al. 1988 and Ol- 
szewski & Aaronson 1985) for this galaxy. The reported in- 
ternal metallicity dispersion in this galaxy is also very low, at 
< 0.2 dex, which introduces little uncertainty in our results. 
Applying our method using the central values of the obser- 
vational parameters we obtain the solid line shown in the 



right panel of Figure (3). Again, the dotted curves represent 
an envelope to a large number of reconstructions obtained 
by changing the observational parameters within their er- 
ror ranges. Of this galaxy we can say that most of its stars 
are older than 12 Gyr. Given the observational errors present 
and the large age of the population of this galaxy, we can not 
draw any inferences on the time structure of the SFR{t), as 
this is totally lost in the noise. The duration of the star form- 
ing episodes can only be concluded to have been < 3Gyr. 
Normalizing through the total luminosity of this galaxy we 
obtain rates of > 400MQ/Myr. The result of applying the 
W test to this galaxy gives 47 ± 5 and 44 ± 4 for the model- 
model and model-data sets, respectively, showing our answer 
to be compatible with the data at better than a 1 cr level. Ol- 
szewski & Aaronson (1985) used a ground based CMD and 
a simple isochrone fitting procedure to conclude the popu- 
lation of this galaxy is uniformly old, with the possibility of 
a 2Gyr spread. 



4.3 Leol 

The observations we obtained for Leol are shown in the 
CMD in the left panel of Figure(4). This contains 11334 
stars, which was reduced to 8691 after exclusion of unsuit- 
able stars (e.g. the red clump), as described for the previ- 
ous two galaxies. The excluded region comprised stars with 
My < 1 and B - > 0.7 and B - 1/ < 0.9, though our 
results are highly insensitive to the details of these cut, as 
the RGB close to the red clump is an age degenerate region. 
The CMD of this galaxy reveals a young MS, but also a RGB 
extending down into the turn off region of a much older MS. 
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Figure 4. 

Left: Observational HR diagram for Leo L Right: Inferred SFR{t) for the central values of the observational parameters, solid line. 
The dashed curves represent the error envelope as defined by the quoted uncertainties in the foreground extinction and distance 
modulus. 



The distribution of stars along the MS region is not uniform, 
and is actually encoding an interesting SFR{t). 

Taking for this galaxy the central values of those given 
by Mateo (1998), [Fe/H] = -L5±0.4,Mv = -11.9, E{B- 
V) = O.OliO.Ol and (m-M)o = 21.99±0.2 (Reid & Mould 
1991, Lee et al. 1993b and Demers et al. 1994), we invert the 
CMD of Leol to obtain the SFR{t) shown by the solid curve 
in the right panel of Figure (4). The dotted curves contain 
all other possible answers compatible with our method and 
the observational parameters taken with their errors. 

The SFR(t) of Leol can be divided into three distinct 
phases, which join continuously with no evidence of a dis- 
crete bursting behaviour. The first of these phases lasted 
from 15 — lOGyr ago, and proceeded at a rate of around 
30MQ/Myr. The following two phases were extended peaks 
of star formation activity centered on ages of 8 and iGyr, 
and having durations and maximum amplitudes of around 3 
and iGyr, at 100 and 150 Mcu/Myr respectively, as shown 
in Figure (4). Any total cessation of the star forming ac- 
tivity can be excluded for ages between 1 and 10 Gyr. As 
time resolution is lost beyond this age, the population be- 
yond 10 Gyr could in principle be a single burst, and appear 
extended because of the observational errors. 

The ground based study of Lee et al. (1993) reaching 
only the youngest turn off points, and subsequent analysis 
of this data set by Caputo et al. (1995) and Caputo et al. 
(1996) using isochrone matching techniques and luminosity 
function methods developed for single age globular clusters, 
revealed the presence of stars of ages 1-3 Gyr. Using more 
recent HST data and comparing to modeled GMDs Gallart 
et al. (1998) describe the star formation history of Leol as 
coming mostly from an episode lasting from 6-2 Gyr ago. 



with the addition of an older component of duration 2-3 
Gyr, in good agreement with our inferred SFR{t) for this 
galaxy. The population box of this galaxy given by Mateo 
(1998) is consistent with our results. 

In this case, the W test gives results showing that our 
inferred SFR{t) is incompatible with the data at a two a 
level. As our inversion method has been extensively tested 
using synthetic CMD's, this result shows the data to be in 
conflict with our input assumptions. This is perhaps not sur- 
prising as this system has a much larger internal metallicity 
spread (0.3 ± 0.1) than the two reviewed previously, which 
together with the errors in the metallicity determination al- 
low for quite a large (~ 1) internal spread. We have thus 
solved for the best fitting single metallicity solution, and 
discovered that internal metallicity dispersions are impor- 
tant. This metallicity spread introduces a time uncertainty 
going from 1-3 Gyr, for ages going from 1-13 Gyr. Solving 
simultaneously for the enrichment and star formation his- 
tories is a problem we shall treat later, as an extension of 
the variational calculus approach. A further possible source 
of the disagreement found between the synthetic CMD re- 
productions of our recovered SFR{t) and the data is the 
difficulty of modeling precisely the error structure present 
in real data, as pointed out by Aparicio & Gallart (1995). 
Although the inversion method itself is highly robust to the 
details of the assumed error structure, the very careful com- 
parison of the W test would pick up any discrepancy be- 
tween the assumed error structure used in generating the 
synthetic CMD's, and that actually present in the data. Fi- 
nally, the presently available isochrones do not take into 
account the relative overabundance of a elements at low 
metallicities, which at some level introduces a slight mis- 
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match between the observed stars and the assumed mod- 
eling. Unfortunately, we can not distinguish between these 
possibilities easily. 

4.4 LeoII 

For this galaxy we obtained 7625 stars, of which we used 
4492 after removing from the analysis the red clump (stars 
with both Mv < 1 and B — V < 0.9), the lower regions 
and some blue stragglers blue wards of V — I = 0.2, which 
do not correspond to any of the evolutionaxy phases in- 
cluded in our isochrones. The full observational CMD is 
shown in the left panel of Figure (5). For the central val- 
ues of the observational parameters of this galaxy as sum- 
marized by Mateo (1998) of [Fe/H] = -1.9 ±0.1, My = 
-8.9, E{B -V) = 0.02 ± 0.01 and (to - M)o = 21.63 ± 0.09 
(Mighcl & Rich 1996, Dcmcrs & Irwin 1993 and Lee 1995) we 
obtain a divergent solution, indicating that the isochrones 
used do not correspond to the stars being analyzed. Chang- 
ing the metallicity to [Fe/H] = —1.75, just marginally out- 
side the errors reported by Mateo (1998), gives a stable con- 
vergence of the method and a significant result. The recent 
study by Mighell and Rich (1996) determined a metallicity 
of [Fe/H] = —1.6 ± 0.25 for this galaxy, consistent with 
what was used here. Our inferred SFR{t) is shown by the 
solid curve in the right panel of Figure (5). Again, the dot- 
ted curves represent an envelope to all alternative recon- 
structions obtained by varying the observational parame- 
ters within their errors. As with Leol, this galaxy shows a 
large internal metallicity spread, which is probably what the 
very sensitive W test detects, also giving an incompatible re- 
sult between the model-model and data-model comparisons 
at more than a two a level. The discussion of this point 
given for the previous galaxy applies also to LeoII, metallic- 
ity spreads will have to be considered for a more accurate 
rendering of the star formation history of this galaxy. 

In this case, we see a gradually rising SFR from 12 Gyr 
to a peak of 160 M© /Myr at 8 Gyr, followed by a somewhat 
more abrupt descent, with star formation activity ending by 
around 6 Gyr ago, as shown by Figure (5) . This result would 
be affected by the internal metallicity spread of 0.3 dex of 
LeoII (Mateo 1998), producing a broadening of around 1-3 
Gyr, see Table (1). Comparing with the study of Mighell and 
Rich (1996), who analyze an HST CMD of LeoII by fitting 
a "fiducial sequence" to the CMD and then comparing it 
to theoretical isochrones to solve for the age of the galaxy 
treated as a single parameter, we find no inconsistencies. 
They obtain an age of 9±1 Gyr for LeoII, with an age spread 
of around 4 Gyr, which is compatible with our results. They 
also report some degree of star formation at ages > lOGyr, of 
which we see no evidence. This discrepancy is probably the 
result of the different methods used in the analysis. Given 
the high age resolution of our method, it is not only the 
median age and a representative value for the spread that we 
obtain, also the shape of the burst is recovered, ruling out for 
example a rectangular burst for this galaxy. Not just the age 
and duration of star formation episodes in these galaxies, but 
also the time structure of them can now be reliably inferred, 
and used in aiding theoretical interpretations of the origin 
of these systems. 



4.5 Draco 

For this Icist galaxy our observational CMD contains 3091 
staxs, most of which are located in the lower, age-degenerate 
region of the diagram, and were thus excluded from the 
analysis, leaving 1210 stars after removing also the hori- 
zontal branch and blue stragglers (stars with both My < 
and B -V < 0.7). Our full observational CMD for Draco 
is shown in the left panel of Figure (6). Mateo (1998) 
summarizes the metallicity and observational properties as 
[Fe/H] = -2.0 ± 0.15, My = -8.8, E{B - V) = 0.03 ± 0.01 
and (to-M)o = 19.58±0.15 (Carney & Seitzer 1986, Lehn- 
ert et al. 1992, Nemec 1985 and Grillmair et al. 1998). Using 
the central values we produced our observational CMD and 
inverted it to obtain the SFR{t) plotted as the solid curve 
in the right panel of Figure (6) . The dotted curves contain 
all variations obtained by shifting our observations within 
the error ranges of the values given by Mateo (1998) for 
E{V- I) and {m- M)v 

Our result for Draco is very similar to what we obtained 
for LeoII, only shifted by 2 Gyr towards younger ages, giving 
a median age of around 7 Gyr. The time structure of the 
SFRit) also differs slightly from that of LeoII in that the 
peak is broader in Draco, having a plateau lasting around 
2 Gyr, rather than a narrow maximum. The presence of a 
low level old component extending beyond 10 Gyr is also 
inferred, although the precise time structure in this region 
is not well restricted by the method. The values given by 
the normalization through the total luminosity are in the 
range of our other galaxies, with the maximum rate being 
of around 110 M^/Myr. Using Saha's W test in this case 
gives 208 ± 9 for the model-model comparison, and 201 ± 8 
for the data-model comparison, which gives us confidence in 
our results, as it shows our SFR{t) is compatible with the 
data at better than a 1 cr level. 

Comparing this result with the recent study of Grillmair 
et al. (1998) we find the two results to be only marginally 
consistent, as they report an age of 10-12 Gyr (for the IMF 
we assumed) for the bulk of the stellax population of Draco 
±2.5Gyr, which they identify as essentially a single age 
event. We note several differences between their approach 
and ours, any of which on its own could bring the two 
results into closer agreement. They use HST data to con- 
struct their observational CMD, which is very similar to 
the one we obtain, no differences are evident at this level. 
Their analysis however differs markedly. They fit a fidu- 
cial sequence to the CMD diagram assumed to be repre- 
sentative of the bulk of the population, adjusting the MS 
and lower RGB regions through an inclusion envelope cri- 
terion, and the bright RGB by eye. This fiducial sequence 
is then compared to theoretical isochrones (VandenBerg & 
Bell 1985) through a maximum likelihood analysis designed 
to find the age of the system, treated as a one parameter 
problem. We note that comparing any fiducial sequence to 
theoretical isochrones will only be a meaningful statistical 
procedure in cases where the underlying SFR{t) is indeed 
a single epoch burst, e.g. in the case of a globular cluster. 
Further, defining any such sequence so that it is a valid sta- 
tistical representation of the underlying SFR{t) is a prob- 
lem that has not been treated yet. Grillmair et al. (1998) 
also noted their isochrones showed systematic inconsisten- 
cies when compared to the stars they were dating, which 
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Figure 5. 

Left: Observational HR diagram for Leo IL Right: Inferred SFR{t) for the central values of the observational parameters, solid line. 
The dashed curves represent the error envelope as defined by the quoted uncertainties in the foreground extinction and distance 
modulus. 



-2 I — . — . — . — . , — . — . — . — . — , — . — . . — . — , — . — . — . — . — , — . — . — I 150 



' 




Figure 6. 

Left: Observational HR diagram for Draco. Right: Inferred SFR{t) for the central values of the observational parameters, solid line. 
The dashed curves represent the error envelope as defined by the quoted uncertainties in the foreground extinction and distance 
modulus. 
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also casts some doubt on their results, as they remarked. 
They also had the difficulty of requiring multiple conver- 
sions between their observational bands and those available 
from theoretical stellar models. 

Our results for this galaxy arc weakened by the as- 
sumption of a single metallicity for the entire population. 
Although this can not be rigorously correct, it essentially 
holds for the previous four galaxies, which show small in- 
ternal metallicity spreads. Draco however, has an internal 
spread of 0.5 dcx (Mateo 1998), which could alter our results 
for this galaxy at the level of 1-3 Gyr. Another possible ex- 
planation to the difference between our results and those of 
Grillmair et al. (1998) is that wo xisod the Padova isochronos 
(Fagotto ct al. 1994, Girardi et al. 1996) rather than those of 
VandenBerg & Bell (1985). Finally, we note that Carney & 
Seitzer (1986) sampled a much larger region of Draco using 
ground based CCD data, and detected multiple turnofFs in 
this galaxy, corresponding to ages between 8 and 15 Gyr. 



5 SUMMARY 

We have used a homogeneous set of observational colour 
magnitude data to study the star formation histories of a 
sample of 5 dSph galaxies, through a non-parametric vari- 
ational calculus maximum likelihood method. We then per- 
formed a detailed statistical analysis to check the accuracy 
of our results for each galaxy, obtaining good results for 
three of our galaxies (Carina, Ursa Minor and Draco), and 
evidence of a systematic difference between our data and 
results for Leol and LeoII, probably due to internal metal- 
licity spreads. We can now compare the results we obtained 
for the different galaxies, with the added consideration of a 
possible extra 1-3 Gyr error margin in the results for Leol 
and LeoII. 

Ursa Minor appears to bo the only essentially "Popula- 
tion 11" system, being characterized by a uniformly old star 
formation history. LeoII and Draco are systems which show 
similar star formation histories, being basically character- 
ized by a single major episode. This lasted in both cases 
around 4 Gyr centered at 8 Gyr, although Draco shows a 
low level extension into much older ages. Leol shows the 
most complex SFR{t), having a small old component of age 
> 10 Gyr, and two later episodes centered at 8 and 3.5 Gyr, 
although the star formation activity did not stop altogether 
between them. It is interesting that the second episode in 
Leol coincides in age with the period of star formation in 
Draco and LeoII. Finally, star formation in Carina highly 
resembles that in Leol in the relative amplitude, duration 
and locations of the two main components. 

Since we used the total luminosities of these galaxies to 
normalize the inferred SFR{t) in physical units, we can de- 
rive other quantities of interest, for example the supernova 
rates as a function of time. The SN type II rates are obtained 
by scaling the total SFR{t) by a factor given by the fraction 
of stars more massive than 8Mq , which for our assumed IMF 
translates 100M©/Mj/r « lSNII/2Myr. The only galaxy 
showing rates greater than I5OM0 /My is Ursa Minor, and 
it is also the only one with a SFR{t) consistent with a sin- 
gle epoch burst. The other four systems, showing extended 
SFR{tys, have rates of always less than 150Mo /Myr. This 
last fact might indicate the presence of a threshold in these 



systems, above which energy input from massive stars into 
the gas component is sufficient to totally disrupt the galaxies 
interstellar medium and end star formation. No character- 
istic timescales of ~ 1 Gyr are evident from the recovered 
SFR(t)'s of these galaxies, with the possible exception of 
Ursa Minor, suggesting that SN type I are not determinant 
in driving the star formation processes in these systems. 

Figure (7) presents the luminosity weighted sum of the 
SFR(t) for all the galaxies we analyzed. This shows the 
star formation activity in the set of these low metallicity 
dSph's to have ended by around 2 Gyr ago, and having been 
relatively steady during the period 3-9 Gyr. For older ages 
we see the average SFR{t) to be essentially dominated by 
the old age Ursa Minor galaxy. Our results as summarized 
by Figure (7) support the calculations presented by Unavane 
et al. (1996) in that the average metal-poor dSph star is of 
intermediate age, and not as old as a "Population 11" halo 
star. A more complete sample would include the Sagittarius 
dwarf, with a mean age of « lOGyr, as well as the much 
larger Magellanic clouds, having ages of < 3Gyr. It seems 
reasonable to suppose the total star formation history for 
the satellites of the Milky Way to show no preferred epoch 
of star formation, as suggested by Tolstoy (1998). 

From these comparisons it is clear that the dSph galax- 
ies of the Milky Way do not form a simple system, and 
straight forward correlations between SFR{t) and other 
present day parameters such as instantaneous galactocentric 
distance or metallicity are not evident, and perhaps not oven 
meaningful. A physical understanding of these systems will 
probably have to consider the complex interactions of these 
systems with the halo of the Galaxy (tidal forces, evolving 
gaseous component, orbital structure etc). It could well be 
the case that the present day sample of survivors actually ex- 
perienced very distinct origins and evolutionary histories (as 
direct studies of their SFR{t) appear to indicate) with little 
in common other than having boon shaped under the domi- 
nating influence of the Milky Way. A ram pressure stripped 
dwarf irregular and a more recently tidally torn fragment of 
the Magellanic clouds could both end up as dSph systems 
today. 

Much more information is needed on these systems bo- 
fore their full potential as tracers of the build up and evo- 
lution of the Milky Way can be realized, once their indi- 
vidual evolutionary histories axe better understood. Stud- 
ies aimed at recovering the full orbital structure of these 
galaxies, through proper motion measurements and poten- 
tial theory reconstructions will yield crucial independent in- 
formation on the evolution and formation of these systems. 
A more complete sampling than the one wo have conducted 
hero is needed to visualize local dSph's fully. Any advance 
in the modeling of advanced stellar phases will also improve 
the use of CMD as tools in galactic evolution studies, as it 
would in principle eliminate the need to remove parts of the 
CMD from consideration. 

Our present sample includes the 5 dSph galaxies hav- 
ing the lowest internal metallicity spreads, and therefore the 
ones for which our present method applies best. Obtaining 
a larger sample and analyzing it through a fully consistent 
and noil parametric statistical method, will require the si- 
multaneous recovering of the enrichment history and the 
SFR{t). Development of such a method will be the subject 
of a future work. We emphasize the need of a homogeneous 
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Figure 7. Sum of the SFR(t)'s of the five galaxies we studied 

sample at all levels of the analysis, together with a fully con- 
sistent statistical inversion method which does not assume 
any a priori structure for the SFR(t) one is solving for, in 
comparative studies of star formation histories. 
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APPENDIX A: THE IMAGES 

The images were recovered from the HST data archive. The 
image numbers, filters and observation dates are given in 
table Al. Table A2 gives the A-to-D gains, exposure times, 
and numbers of images for the dSph fields. 
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Al Retrieving the data and image preparation 

For each of the five dSph's, there wore a number, n, of images 
taken, in both a long band (F814W, corresponding closely 
to Johnson's I), and a short band (F555W or F606W, cor- 
responding to Johnson's V). 

For each dSph, there are a sot of long exposure 'V' and 
T data files. Each data file contains the output from 4 detec- 
tors - the planetary camera (PC) , and the three Wide Field 
cameras (WF2, WF3 and WF4). Each of these detector im- 
ages is 800x800 pixels in size, of which typically 730x730 
are usable. 

Data treatment was carried out in the IRAF environ- 
ment. Image combination was carried out using the STSDAS 
package, and photometry using the DAOPHOT package. 

In order to remove the severe cosmic ray effects in each 
image, the n images were combined by taking the mean value 
for each pixel position, after the rejection of values either too 
high or too low with respect to local variations. The task 
used here was 'crrej' in the STSDAS package. Each image 
was also cropped to a 730x730 image by rejecting the first 
60 rows/columns and the last 10. 

In the following analyses, only the data from the three 
WF cameras were used. The pixel scale in these WF images 
is 0.10" per pixel. 

A2 Source extraction 

Sources were extracted to 2a above the mean background. 
Point spread function (PSF) fitting photometry was carried 
out by selecting isolated stars to define a PSF. Sources with 
a fit at > 1-5 were rejected, and the magnitudes used 
were aperture magnitudes with a radius of 2 pixels. (0.20"). 

A3 Galactic contamination 

At faint magnitudes in optical wavebands, external galax- 
ies can constitute a major contaminant in number counts. 
Williams et al. (1996), based on HDF (Hubble Deep Field) 
galaxy counts find ~2xl0'' galaxies per square degree 
brighter than V=26. For the area of the three WF detectors, 
this corresponds to, on average, 250 contaminating galaxies 
in the field of view. Compared with the many thousands of 
stars in the images, this contamination is small. 

Nevertheless, the high resolution of the HST allows 
the separation of galaxies from stars more reliably than 
for ground based work, where image resolution is necessar- 
ily lower and the distinguishing stars from galaxies is less 
straight forward . 

The cut in which we use eliminates the majority of 
galaxies (and remaining cosmic ray events) because these 
will in general be fitted poorly by the PSF. 

A4 Magnitude corrections 
A4.1 Aperture Corrections 

An aperture correction to render these magnitudes equiva- 
lent to those which would be obtained by use of a 0.5" radius 
aperture (see Holtzman et al. 1995) was found by selecting 
bright, unsaturated stars, in each of the 'V and T bands, 
and for each detector and each dSph separately. The mean 



differences between the magnitude using a 5.02 pixel diam- 
eter aperture (0.5" radius) and a 2 pixel diameter aperture 
were used to correct all magnitudes. 

A4-2 A to D gain correction 

All the observations for the dSph's considered here were 
taken through bay 4 (see Holtzman et al. 1995), which means 
that the Analogue-to-Digital gain is only 7.0, rather than 
the standard value of 14.0. Due to some unshared electron- 
ics, this necessitates a different correction for each of the 
WF fields, as indicated below: 

Wide Field Am 

2 0.754 

3 0.756 

4 0.728 



A4-3 Geometric correction 

The WF cameras have geometric distortions which arise 
mainly from elements in the optical path. As a consequence, 
the effective pixel areas, in square angular measure, vary 
systematically across each WF detector. We make a param- 
eterization of the data from the figures of Holtzman et al. 
(1995), and apply that as a correction. The correction is well 
represented by a quadratic function of distance from the cen- 
tre of the detector, and never exceeds 0.04 magnitudes at the 
edge of the detectors. 

We use Am = -1.897 x 10^^ r- 1.208 x lO"'^ r^ where 
r is the distance in pixels from the position (400,400) on the 
detector. 

A4.4 Charge Transfer efficiency correction 

The readout of the CCD detectors requires the transfer of 
charge through successive rows of the detector. As a con- 
sequence, the signal from the last rows to be read are di- 
minished because of the loss of charge during transfers. The 
correction is a maximum of 0.04 magnitudes at the final row. 
We use Am = — 0.04(j//800) where y is the row number. 

A4-5 Corrections for Leo I 

WFPC2 data taken before 23rd April 1994 was taken at 
a detector temperature of — 76°C rather than — 88°C. The 
change reduced CTE effects, and IR zeropoint problems. 
The Leo I observations were made before this change, and 
we use a linear ramp of size 0.12 magnitudes to correct for 
the CTE effects as recommended by Holtzman et al. (1995), 
and additionally, an additive offset of 0.05 magnitudes in I 
to correct for the higher zero point. 

A5 Conversion to standard V and I 

Reddening must be taken into account before corrections are 
made to V and I. We use the published values of reddening, 
which are fairly small for all these galaxies {E[B — V) < 
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Table A3. The synthetic transformations from F555W and 
F606W to V, and from F814W to I, taken from Holtzman et 
al. (1995) 



Filter 


Constraint 


Output Band 


ao 


ai 


0-2 


F555W 




V 


21.729 


-0.051 


0.009 


F606W 


(V-I)<1.0 


V 


22.093 


0.254 


0.012 


F606W 


(V-I)>1.0 


V 


22.883 


-0.247 


0.065 


F814W 


(V-I)<1.0 


I 


20.838 


-0.012 


-0.006 


F814W 


(V-I)>1.0 


I 


20.920 


0.028 


-0.124 



0.08), to make corrections to the magnitudes before applying 
the transformations given below. 

We use the extinction values tabulated by Holtzman et 
al. (1995) for the filter F555W and F814W, and an estimate 
based on these for the F606W filter: 

Af555w = 3.026£(B - V) 

Af814w = l.S25E{B-V) 

Af606W = 2.75E{B - V) (estimate) 

We subsequently apply synthetic transformations given 
by Holtzman et al. (1995), of the form 

Output Band — rriraw + ao + ai(V— I) + a2(V— I)^ 



where the coefficients are given in table A3. 

rriraw is the output aperture magnitude using an aper- 
ture of 2 pixels radius, corrected as indicated above to an 
aperture of radius 0.5". 

The above formulae were iterated until no further sig- 
nificant change in V or I occurred. 

Finally, the reddenings which had been removed are re- 
stored using: 

Av = 3.10E{B ~V) 

and 

Ai = 1.83E{B-V) 

The reddenings used are refined iteratively after the fit- 
ting of isochrones to the MS region. Note, however, that 
because the reddenings are small, and the reddenings are 
restored to the data afterwards, even if the initial guess for 
the reddening is wrong by a substantial amount, very little 
effect is seen in the final photometry, (e.g. a change of the 
reddening by 0.1™ changes photometry by less than 0.005™.) 



A6 Systematics 

As Holtzman et al. (1995) point out, there remain some as- 
pects of the photometric calibration of WFPC2 which are 
uncertain. For example, Holtzman et al.(1995) note discrep- 
ancies of ~ 0.05 magnitudes between long and short expo- 
sures, which are not understood. Several more minor sys- 
tematic and random effects at the level of a few percent 
(corresponding to a few hundredths of a magnitude) are not 
well understood. Furthermore, the conversion that has been 
made here to standard V and I colours, introduces more 
minor uncertainties. It must be noted that the systematic 
effects in zeropoints may be as large as 0.1 magnitudes. Any 
such uncertainties appear as an ofTset primarily in adopted 
distance modulus, and would affect most CMDs in the same 
way. Differentially, any effects should be minimal. 



